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ABSTRACT 


Continuous- time  models  of  oscillator  phase  noise  x(t)  usually  have 
stationary  nth  differences,  for  some  n.  The  covariance  structure  of  such  a 
model  can  be  characterized  in  the  time  domain  by  the  structure  function: 
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Although  formulas  for  the  special  case  D2(0;t,t)  (the  Allan  variance  times 
2x2)  exist  for  power-law  spectral  models,  certain  estimation  problems  require 
a  more  complete  knowledge  of  (0) . 

We  exhibit  a  much  simpler  function  of  one  time  variable,  D(t),  from 
which  (0)  can  easily  be  obtained  from  the  spectral  density  by  uncomplicated 
integrations.  Believing  that  D(t)  is  the  simplest  function  of  time  that 
holds  the  same  information  as  (0),  we  call  D(t)  the  fundamental  structure 
function. 

We  compute  D(t)  for  several  power-law  spectral  models.  Two  examples  are 
D(t)  =  K j  t i 3  for  random  walk  FM,  D(t)  -  Kt2  ln|t|  for  flicker  FM.  Then,  to 
demonstrate  its  use,  we  exhibit  a  BASIC  program  that  computes  means  and  vari¬ 
ances  of  two  Allan  variance  estimators,  one  of  which  incorporates  a  method  of 
frequency  drift  estimation  and  removal.  Except  for  a  one-line  function  defi¬ 
nition  of  D (t) ,  the  program  is  independent  of  the  phase  noise  spectrum.  The 
outputs  were  used  for  assigning  confidence  intervals  to  the  results  of  recent 
hydrogen  maser  performance  tests  at  JPL. 

THE  STRUCTURE  FUNCTION 

The  purpose  of  this  paper  is  to  demonstrate  an  easy  way  of  computing  a 
class  of  time-domain  oscillator  stability  measures  known  collectively  as  the 
structure  function.  Let  the  phase  of  an  oscillator  with  nominal  frequency  vq 
be  modeled  by 


2'ttVq  ( t  +  x  ( t ) )  , 

where  x(t)  is  a  random  process  representing  the  "phase  time"  of  the 


*The  research  described  in  this  paper  was  carried  out  by  the  Jet  Propulsion 
Laboratory,  California  Institute  of  Technology,  under  contract  with  the 
National  Aeronautics  and  Space  Administration. 
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oscillator.  The  most  restricted  version  of  the  structure  function,  discussed 
by  Lindsey  and  Chie  [1],  is 

Dn(T)  =  E[A*  x(s)J2  ,  (1) 

where  AT  is  the  backwards  difference  operator  defined  by  ATf (t)  =  f (t)  -  f (t-r) , 
and  E  denotes  the  mathematical  expectation.  It  is  to  be  understood  that  the 
nth  difference  process  A1^  x(t)  is  wide-sense  stationary,  so  that  (1)  does  not 
depend  on  s.  For  n  =  2  we  obtain  the  Allan  variance  as 


o2(t) 

y 


D2(t)  . 


A  more  general  version  of  the  structure  function  is 

Dn(t;  =  E  i"  x(s+t)  A^  x(s)  ,  (2) 

1  2 

which  was  used  by  Yaglom  [2]  as  a  basis  for  a  theory  of  processes  with 
stationary  nth  differences.  For  the  application  to  be  given  at  the  end  of 
this  paper,  we  shall  need  the  yet  more  general  version 

D(t;  a.,b)  =  E  Aax(s+t)  A^x(s)  ,  (3) 

where 


a  -  (a1,...,an),  Aa  -  Aa  . ,.Aa  , 

—  In 

and  likewise  for  _b,  A^. 

By  using  the  transfer  function  1  -  e  of  the  operator  AT}  one  can 
express  all  of  these  quantities  as  integrals  involving  the  two-sided  spectral 
density  Sx(u>)  of  the  process  x(t) ,  which  is  assumed  to  have  stationary  nth 
differences  for  some  n.  For  example,  if  the  long-term  frequency  drift  rate  is 
zero ,  then 


D2e) -/ n  -  e-^|  v.)£ 


CO 

=  16  /  sin  Sx(w)  ~  , 

—  CO 

D(t;  a,b,c,d)  =  E  A^A^  x(s+t)  AcAd  x(s) 

oo 

fa-  e-ia)a)(l  -  e-iu)b)(l  -  ela,c)(l  -  eiud)  elwt  Sx«o) 


(4) 


(5) 


(This  is  (3)  with  n  =  2.) 
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Because  of  the  importance  of  the  Allan  variance  as  a  stability  measure, 
the  integral  (4)  has  been  evaluated  in  closed  form  for  power-law  spectral 
models 

S  («o)  -  K  |  m  | a_2,  a  =  -2,  -1,  0,  1,  2, 

with  a  high-frequency  cutoff  when  ail  ("PM"  noises)  [3].  The  expressions 
(5)  and  (3)  are  more  complex,  yet  we  shall  show  how  to  compute  any  of  these  in 
two  easy  steps:  1)  By  means  of  simple  integrations  involving  only  the  spec¬ 
tral  density,  evaluate  a  single  function  D(t)  of  one  time  variable;  2)  apply 
a  certain  difference  operator  of  order  2n  to  D(t).  This  function  D(t), 
depending  only  on  the  spectral  density,  will  be  called  the  fundamental  struc¬ 
ture  function  because,  merely  by  taking  linear  combinations  of  values  of  D(t), 
one  can  generate  the  entire  covariance  structure  of  the  differences  of  x(t). 

STATIONARY  PHASE 

If  x(t)  is  stationary,  then  D(t)  turns  out  to  be  just  the  autocovariance 
function : 

00 

D(t)  =  Cov(x(s+t)  ,x(s))  =  f  ela)t  S  (to)  4^  .  (6) 

X  Z7T 

—oo 

Indeed,  by  expanding  the  differences  in  the  middle  expression  of  (5),  one 
can  show  directly  that 

D(t;  a,b,c,d)  =  Aq  A_c  A_d  D(t)  .  (7) 

It  is  just  this  form  of  expression  that  will  be  extended  to  the  nonstationary 

situation,  where  S  (m)  is  not  integrable  near  zero  frequency. 

X 

THE  FUNDAMENTAL  STRUCTURE  FUNCTION  IN  GENERAL 

Let  us  be  given  a  process  x(t)  with  stationary  nth  differences.  It  is 
known  [2]  that  x  has  a  spectral  density  sx(w)  with  the  properties 


< 

/ 


S  (to)  dco  < 
x 


S  cok  S  (oo)  doo  < 
0 


oo 


(8) 


for  some  integer  k,  0-k  -  2n.  (We  ignore  the  more  general  case  of  a  spectral 
distribution  function.)  Here  is  one  way  to  compute  the  fundamental  structure 
function  D(t) :  Define  the  function 

BU>  »  /V“Z  (iU)kSx(U>  ff  (9) 


283 


on  the  upper  half-plane  Im  z  >  0.  Let  C(z)  be  any  kth  integral  of  B(z),  i.e., 
any  function  such  that  c(k)  (z)  =  B(z)  on  Im  z  >  0.  It  will  always  be  found 
that  C(z)  can  be  extended  to  the  real  line  by  continuity.  Then  let 

D(t)  =  2  Re  C(t) ,  t  real  (10) 

If  k  >  0,  then  D(t)  is  not  unique;  any  polynomial  of  degree  <  k  may  be  added 
to  it.  If  x  is  stationary  then  k  =  0,  and  D(t)  reduces  to  the  covariance 
function  (6). 

RANDOM  WALK  FM 

To  see  how  easy  it  is  to  carry  out  the  above  procedure  for  power-law 
noise  types,  consider  the  spectral  density  S  (w)  =  and  take  k  =  4.  Then 

X 

oo 

\  C  i“z  dm  1  nx 

B(Z)  =JQe  2?=  "  2^17  (Im  z  >  0)  * 

a  fourth  integral  of  which  is 

C(z)  =  -  -^-r  — — (Im  z  >  0)  . 

2tti  6 

(One  can  throw  out  a  polynomial  of  degree  less  than  4.)  Taking  In  z  = 
ln[z|  +  i  Arg  z  for  Ii  z  >  0  gives 

D(t)  =  2  Re  C(t)  =0  (t  >  0) 

=  -t3/ 6  (t  <  0). 

3 

Adding  t  /12  to  this  gives  the  alternate  form 

D(t)  =  1 1 1  3/12  . 


OTHER  POWER-LAW  NOISE  TYPES 

We  now  give  D(t)  for  the  power-law  spectrum 


S  (m)  =  K  u 
x  a. 


a- 2 


h 

a 

2(2TT)a 


The  forms  for  the  constant  Ka  and  the  power  of  m  link  the  results  to  an 
accepted  notation  for  the  noise  spectrum,  namely 


Sy(f) 


h  f 
a 


a 


[3],  where  y  =  dx/dt,  f  =  co/(2ir),  and  Sy(f)  is  the  one-sided  spectral  density 
of  y.  Fractional  values  of  a  are  allowed.  We  can  take  k  to  be  the  integer 
part  of  2  -  a  (refer  to  (8)). 
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Case  1. 


a  <  1,  not  an  odd  integer. 


D(t)  "  Ka  2r(2-a)cos(na/2)  '  (11) 

This  case  includes  white  FM  (a  =  0)  and  random  walk  FM  (a  =  -2). 

Case  2 *  ot  —  —  1 ,  — 3 ,  —  5 ,  •  *  . 


il-a 


D(t)  ■  (_1)(3-0/2  . 


(12) 


This  case,  which  includes  flicker  FM  (a  =  -1),  is  actually  the  easiest  of  all. 
Case  3.  a  =  1  (Flicker  PM). 

Here  we  use  the  spectrum 

Sx(oj)  =  K-^lwl-1  exp  (-jto/aJh|)  . 


The  exponential  high-frequency  cutoff  yields  the  elementary  result 

D(t)  =  -  ^  ln(t2  +  1/rn2)  ,  (13) 

whereas  the  usual  rectangular  cutoff  yields  a  cosine  integral. 

Case  4.  a  >  1  . 

This  phase  noise  process  (provided  with  a  high-frequency  cutoff)  is 
stationary;  thus  D(t)  is  just  the  autocovariance  function  of  x. 

DERIVING  THE  STRUCTURE  FUNCTION  FROM  D(t) 

Having  the  fundamental  structure  function  D(t)  in  hand,  we  now  show  how 
to  use  it.  Let  x(t)  be  a  noise  with  spectral  density  Sx(u)) ,  k  the  smallest 
integer  satisfying  (8) ,  and  tig  the  smallest  integer  such  that  2ng  ^  k.  The 
differences  of  x  of  order  ng  and  higher  are  stationary,  and  to  simplify 
matters  let  us  assume  that  they  are  ergodic  and  have  mean  zero.  For  ng  =  2, 
this  means  that  an  oscillator  with  phase  time  x(t)  has  no  long-term  average 
frequency  drift.  Our  main  result  is  the  following  formula  for  the  general 
structure  function  (3)  of  order  n  ^  ng: 

D(t;  a,b)  =  A  A  D(t) .  (14) 

3-  l) 

The  proof,  which  has  appeared  elsewhere  [4],  starts  from  the  spectral 
integral 
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(15) 


D(t; 


A  A  . 

a  -b 


ellot  s  <„) 

x  2tt 


In  which  the  differences  operate  on  eiwt  as  a  function  of' t.  (Refer  to  (5) 
for  the  case  n  *  2.)  One  cannot  simply  pull  the  difference  operators  outside 
the  integral,  for  the  Fourier  transform  integral  of  Sx(w)  diverges  (unless 
k  =  0)  if  the  singularity  at  w  =  0  is  not  neutralized.  One  way  to  do  this  is 
through  the  function  B(z)  given  in  (9).  Another  embodiment  of  D(t)  is  given 
directly  by  the  formula 


D(t) 


GO 

[cos  cot  - 


1  +  CO 


2n 


n-1 

s 


J“0 


(iC0t)2Ji 
(2j) !  J 


Sx<“>  f 


(16) 


where  n  n^  . 


MOMENTS  OF  STABILITY  ESTIMATORS 

The  remainder  of  this  paper  presents  a  nontrivial  application  of  the 
preceding  theory.  Let  the  phase  time  x(t)  have  a  long-term  frequency  drift 
component : 

x(t)  =  Xq ( t )  +  -|ct2  , 

where  c  is  the  constant  rate  of  frequency  drift,  and  xg(t)  is  a  Gaussian 
process  whose  2nd  differences  are  stationary,  ergodic,  and  have  mean  zero. 
Observing  x(t)  for  0  $  t  -  T,  we  wish  to  estimate  c  and  the  quantities 


c2(x)  =  -i*  E  [A2  x 
y  2t 


(t)]' 


(Gross  Allan  variance)  , 


(t) 


x0(t)]' 


(Net  Allan  variance)  , 


which  are  the  theoretical  Allan  variances  before  and  after  removal  of  drift. 

To  define  and  manipulate  the  estimators,  it  is  necessary  to  set  up  some 
notation  as  follows:  Write 

C(a,b, t)  -  A^xU)  . 

(This  has  nothing  to  do  with  (10),)  Then  EC  =  c  for  any  a,b,t.  A  particular 
one  of  these  is  used  for  estimating  c,  namely 

c  -  C(xc,T-Tc,T)  (17) 
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where  tc  =  T/6.29,  a  value  chosen  to  minimize  the  variance  of  c  in  the 
presence  of  flicker  FM  noise  [5].  The  interpretation  of  (17)  is  that  the 
estimated  drift  rate  equals  the  average  frequency  near  the  end  of  the  record, 
minus  the  average  frequency  near  the  beginning,  divided  by  the  intervening 
time  T  -  Tc. 

Let  us  introduce  the  additional  quantities 

Cj  =  C(T,T,jx),  cT  =  C(t,T-t,T)  ,  (18) 

where  j  is  an  integer.  To  estimate  (2/t2)  cr.|(t),  we  use  the  unbiased 
statistic 


(19) 


where  it  is  now  assumed  that  T/t  is  an  integer  m.  This  is  just  the  usual 
Allan  variance  estimator  involving  the  sum  of  squares  of  adjacent  second 
differences.  (The  unusual  scale  factor  2/t^  is  convenient  for  these  calcu- 
lations.)  To  estimate  (2/tz)  oZq(t),  we  use 


V 


0 


1 

m-1 


(20) 


Since 


m 


we  have 


V0  = 


V  —  2cc  +  c 


(21) 


Our  goal  is  to  compute  the  mean  and  variance  of  Vq,  a  biased  estimator. 
First,  we  see  that  Vq  does  not  depend  on  the  true  drift  rate  c,  and  so  for 
this  purpose  we  can  take  c  =  0.  Then 

EV  =  (2/t2)  o20(t)  ,  (22) 
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Var  V  = 


1 


(23) 


(m-1) 


2  EE  Cov(crck)  * 


j  k 


EVq  =  EV  -  2Ecct  +  Ec  , 


(24) 


Var  VQ  =  Var  V  +  4  Var(cc^)  +  Var 


-  4  Cov(V,cct)  +  2  Cov(V,c2) 
-  4  Cov(cc  ,c2)  . 


(25) 


Because  x(t)  is  Gaussian,  the  fourth  moments  in  (25)  are  determined  from 
the  second  moments  of  c,  c.,  c^,  and  the  general  formula 

Cov(uv,xy)  =  Eux  Evy  +  Euy  Evx  ,  (26) 

where  u,  v,  x,  y  are  zero-mean  random  variables  with  a  joint  Gaussian 
distribution.  In  turn,  the  required  second  moments  are  just  special  cases 
of  the  structure  function  formula  (14),  which  here  takes  the  form 

E  C(a,b,s  +  t)  C(c,d,s)  =  VbA-cA-dD(t)  ’  (27) 

where  D(t)  is  the  fundamental  structure  function  of  the  phase-time  noise 
x(t).  This  expression,  when  written  out,  is  a  sum  of  16  terms  involving  D. 
(See  line  1040  of  the  BASIC  program  in  Fig.  1.) 

A  BASIC  PROGRAM 

The  tedious  but  elementary  task  of  computing  (23)  -  (25)  by  means  of 
(26)  and  (27)  is  carried  out  by  the  BASIC  program  shown  in  Fig.  1.  The  main 
point  to  observe  about  this  program  is  that  the  whole  algorithm  is  independent 
of  the  phase  noise  spectrum,  which  enters  only  through  the  one-line  function 
definition  of  D(t)  in  line  1030  (and  the  print  statement  130).  The  only 
requirement  on  Sx(w)  is  that  the  integer  k  in  (8)  be  at  most  4.  This  guaran¬ 
tees  that  x(t)  has  stationary  second  differences.  Because  of  the  way  the  out¬ 
puts  are  scaled,  a  multiplicative  constant  in  D(t)  can  be  neglected;  thus  we 
can  let  D(t)  =  |t|2  for  random  walk  FM.  For  flicker  FM,  line  1030  becomes 

DEF  DD(T)  =  T*T*LOG(ABS(T)  +  (T=0))  . 

Line  1040  evaluates  (27),  and  the  subsequent  definitions  evaluate  the  required 
special  cases.  Note  that  the  length  of  the  test  run  is  used  as  a  unit  of 
time.  It  must  be  admitted  that  computational  efficiency  has  been  sacrified  to 
gain  ease  of  coding. 

Figure  2  shows  the  output  of  the  program  for  random  walk  FM.  The  column 
MEAN(NET)  gives  the  expected  value  of  (t^/2)Vq  relative  to  the  true  net  Allan 
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variance  c'yo(T) •  The  variances  of  V  and  Vq,  the  estimators  of  gross  and  net 
Allan  variances,  are  presented  in  terms  of  "degrees  of  freedom"  (DF) ,  defined 
for  a  positive  random  variable  X  by 


DF 


2 (EX) 2 
Var  X 


as  if  X  had  a  x2  distribution.  The  DF(GROSS)  result  is  valid  only  if  the  true 
drift  has  been  subtracted  from  the  phase  data.  When  T/t  =  1,  the  gross  and 
net  DF  are  both  1,  since  both  V  and  Vq  are  squares  of  mean-zero  Gaussians. 

CONFIDENCE  INTERVALS  FOR  ALLAN  VARIANCES 


Following  Howe  et  al.  [6],  we  use  the  results  shown  in  Fig.  2  to  assign 
confidence  intervals  based  on  the  x2  distribution,  usually  with  a  fractional 
DF.  The  author  has  not  estimated  the  error  caused  by  pretending  that  V  and  Vq 
are  proportional  to  x2  variables.  Random  walk  FM  was  used  as  a  noise  model  for 
t  >  IQ4  s.  Fig.  3  shows  the  estimates  of  gross.  Allan  variance  (before  drift 
removal)  and  net  Allan  variance  (after  drift  removal)  for  a  72-day  test  run  of 
a  pair  of  hydrogen  masers  at  JPL. 

Different  drift  estimators  are  used  in  Figs.  3a  and  3b.  In  Fig.  3a,  we 
use  a  value  of  drift  measured  by  retuning  the  masers  over  a  period  much  longer 
than  the  72-day  test  run.  Regarding  this  value  as  close  to  the  "true"  value 
c,  and  using  it  in  (21)  in  place  of  c,  we  can  use  the  gross  DF  numbers  from 
Fig.  2  to  assign  confidence  intervals  to  net  Allan  variance.  In  Fig.  3b,  the 
estimated  value  c  was  used,  so  that  one  must  use  the  net  mean  and  DF  numbers 
to  compute  the  confidence  intervals.  The  negative  bias  of  Vq  pushes  the  con¬ 
fidence  intervals  up;  in  fact,  for  T/ x  =  2  (one  sample  of  second  phase  differ¬ 
ence),  the  bias  is  so  great  that  the  90%  confidence  interval  does  not  contain 
the  estimate  itself.  This  is  pushing  things  too  far,  perhaps. 

The  basic  problem  is  that  one  cannot  remove  estimated  drift  (via  (21), 
for  example)  without  also  taking  a  bite  out  of  the  long-term  random  fluctua¬ 
tions.  The  outputs  of  the  above  computer  program  for  oja"^  phase  noise 
(-2.5  4  a  4  0)  support  the  conjecture  that  the  bias  of  the  net  Allan  variance 
estimator  Vg  is  always  negative. 
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100  ! MOMENTS  OF  ALLAN  VARIANCE  ESTIMATORS 
110  GOSUB  1000 
120  MC-6 . 29 

130  PRINT  \ ' NOISE  TYPE  «  RANDOM  WALK  FM# 

140  PRINT  'MEAN  A  D.F.  FOR  GROSS  A  NET  A.V.'W 

'T/TAU' , 'MEAN (NET) ' , 'D .F . (GROSS) ' > 'D.F. (NET) ' 

ISO  FOR  M«2  TO  9,10  TO  18  BY  2,20  TO  SO  BY  S 
160  GOSUB  1080 

170  DF«2/VARV,DF0«2#EV0*EV0/VARV0 
180  PRINT  M,  EV0 ,DF ,DF0 
190  NEXT  M 
200  STOP 

1000  ! DEFINITIONS 

1030  DEF  DD(T)»ABS(T#T#T)  ! FUNDAMENTAL  STRUCTURE  FUNCTION 

IFOR  RANDOM  WALK  FM 

1040  DEF  MM(A,B,C,D,T)  =  (DD(T)  ~DD  ( T~A  >  -DD(  T-B)-*DD  (T+C)~DD(T+D> 

+DD ( T-A-B ) +DD ( T-A+C ) +DD ( T-A+D ) +DD(T-B+C )+DD(T-B+D)+DD(T+C+D) 
-DD  ( T-A-B+C )  -DI)  ( T-A-B+D  )  -DD  ( T-A+C+D) ~DD  ( T-B+C+D  > 
•»-DD(T-A~B+C+D>  )/(A#B*C*D) 

1050  DEF  FF<A,B>=*MM(A,1~A,B,1“B,0> 

1060  DEF  GG<A,B,J>=*MM(A,l-A,B,B,i-J*B) 

1070  DEF  HH(A, J >=MM < A, A , A , A , J#A) 

1073  RETURN 

1080  ! COMPUTATION  SUBROUTINE 

! INPUTS:  M=T/TAU>  MC=T/TAUC  <~6.29> 

IOUTPUTS:  EV0/A.V.,  VAR  ( V ) /A  .  V  .  **2 ,  VAR  ( V0  >/A  .  V  .  *#2 
EV0  VAR  V  VAR  V0 

1090  U-i/M,UOi/MC  !U  STANDS  FOR  TAU 
1100  FUU~FF(U,U>  ,FCUtsFF(UC>U)  ,FCC“FF(UC,UC) 

1110  HU0=-HH(U,Q) , HU02=HUQ#HUQ  !HU0  IS  SCALED  ALLAN  VARIANCE 
1120  ! 

1130  EVO=HU0-2*FCU+FCC 

1140  Ml*=M~i,X»0\  H»HH(U,  J)  ,X*=X+(Mi-J)#H*H  FOR  J»i  TO  M-2\ 
X»2*X+Mi#HU0 2 , VAR V=X*2/ ( M 1 *M 1 ) 

11S0  VARCU>=FCC*FUU+FCU*FCU,  VARCC«2*FCC*FCC 

1160  X«0\  X*X+GG(UC,U,J)*GG<U,U, J>  FOR  J*2  TO  M\  C0VVCU=2*X/Mi 
1170  X»0\  G»GG(UC,U,J)  , X-X-^GHtG  FOR  J=2  TO  M\  C0VVCC*2*X/M1 
1180  CVCUCC=>2*FCC#FCU 
1190  ! 

1200  VAR V0=VARV+4*VARCU+VARCC-4*CQVVCU+2#COVVCC-4*CVCUCC 
1210  RETURN  WHERE  EVO»EV0/HU0 , VARV»VARV/HU02, VARVO=VARV0/HUO2 


Figure  1.  A  BASIC  program  for  computing  mean  and  variance 
of  an  estimator  of  net  Allan  variance. 


NOISE  TYPE  *  RANDOM  WALK  FM 
MEAN  &  D.F.  FOR  OR  OSS  6  NET  A.V. 


T/TAU 

MEAN (NET) 

D.F. (CROSS) 

D.F. (NET) 

2 

.11213718 

1 

1  .  0000011 

3 

,4131003 

1 . 882353 

1 . 2011257 

4 

. 56608639 

2.7692308 

1 . 9797428 

5 

. 6S837896 

3.6571431 

2 . 8213698 

6 

.72007427 

4 . 5454549 

3 . 6927653 

7 

.76417726 

5 . 4339623 

4.5779951 

8 

.7970189 

6.3225806 

5 . 4662905 

9 

.82222714 

7.2112679 

6 . 3534235 

10 

.84209356 

8. 1000005 

7 .2390502 

12 

.87125838 

9.8775517 

9 . 0083684 

14 

.89153524 

11 . 655173 

.10.777728 

16 

.90639572 

13.432836 

12.546251 

18 

. 91772997 

15.210527 

14.314574 

20 

. 92664775 

16.988236 

16.084209 

25 

.9423454 

21 . 432559 

20.511747 

30 

. 95254386 

25 . 876923 

24 . 943548 

35 

.9596919 

30.321313 

29  .378236 

40 

.96497606 

34 .765708 

33.814985 

45 

.96903914 

39.210128 

38.253179 

50 

.97225997 

43.654528 

42.692561 

Figure  2.  Output  of  the  program  of  Fig.  1  for 
random  walk  FM  noise. 
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QUESTIONS  AND  ANSWERS 


None  for  Paper  #12 
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